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• The possibihty to extract relevant information on spectroscopic factors from (e,e'p) reactions at 

I high is studied. Recent ^^0(e,e'p) data at — 0.8 (GeV/c)^ are compared to a theoretical 

, ^ . approach which includes an eikonal description of the final-state interaction of the proton, a micro- 

ti ' scopic nuclear matter calculation of the damping of this proton, and high-quality quasihole wave 

, functions for p-shell nucleons in ^^O. Good agreement with the = 0.8 (GeV/c)^ data is obtained 

■ when spectroscopic factors are employed which are identical to those required to describe earlier low 

' experiments. 
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' During the last fifteen years a quiet revolution has taken place in the perception of the nucleus. During this period 
fSI the study of the (e,e'p) reaction has clarified the limits of validity of the mean-field description of nuclei. In particular, 
' absolute spectroscopic factors have been obtained for the removal of protons from many nuclei [Q-Q . The qualitative 
, features of the strength distribution suggest a considerable mixing between single-hole states and more complicated 
' configurations, like two-hole— one-particle states. The resulting fragmentation pattern of the single-particle (sp) 
I [ strength exhibits a single peak carrying about 65-70% of the strength for states in the immediate vicinity of the 
^ ■ Fermi energy, while more deeply bound orbitals display a strongly fragmented distribution reminiscent of a complex 
(— ( sp energy. In addition to these fragmentation features, an important depletion of the sp strength has been established 
• • 1 which is associated with ground-state correlations induced by strong short-range and tensor correlations This 
' leads to an overall reduction of the sp strength for all mean-field orbits in all nuclei by 10-15% This theoretical 

, result has recently been confirmed also for deeply-bound orbits in an (e,e'p) experiment on ^"^^Pb in a wide domain 
^ of missing energy and momentum |^,^. 

Ci The analysis of the (e,e'p) data has relied on the Distorted Wave Impulse Approximation (DWIA) for both the 

Coulomb distortion of the electron waves (in heavy nuclei) and the outgoing proton ||lO|-p^ . The proton distortion is 
described in terms of an optical potential required to describe elastic proton scattering data at relevant energies 
There is some uncertainty related to this treatment since elastic proton scattering is considered to be a surface reaction 
and no detailed information is obtained related to the interior of the nucleus. This uncertainty gives rise to an 
estimated error of about 10%. Such an estimate may be inferred by considering the difference between the relativistic 
and nonrelativistic treatment of the proton distortion. It is shown in Ref. |lj] that this difference is essentially due 
to the reduction of the interior wave function in the relativistic case. This feature can also be generated by including 
a reasonable amount of nonlocality in the optical potential ]l4| . 



A serious challenge to the interpretation of (e,e'p) experiments was recently published in Refs. 1 15 1^ . This challenge 
consists in questioning the validity of the constancy of the spectroscopic factor as a function of the four-momentum 
(squared), Q^, transferred by the virtual photon to the knocked-out nucleon. In Ref. |l^ a conventional analysis of the 
world's data for (e,e'p) experiments on ^^C at low generated results consistent with previous expectations. Data at 
higher values of were then analyzed within the framework of a theoretical model which employs Skyrme-Hartree- 
Fock bound-state wave functions for the initial proton, a Glauber-type description of the final-state interaction of the 
outgoing proton, and a factorization approximation for the electromagnetic vertex ||l5|,[l6|]. Within the framework of 
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this theoretical description, spectroscopic factors were obtained which increase substantially with increasing for 
the ^^C nucleus. 

The spectroscopic factor is a many-body quantity defined without reference to a probe. For valence hole states in 
nuclei it simply represents the probability for the removal of a nucleon with prescribed quantum numbers from the 
ground state of the target while ending up in a state of the nucleus with one particle less. In the conventional analysis 
of the experimental data these quantum numbers involve the trivial values of parity and total angular momentum but 
also require the corresponding wave function to be a solution of a Woods-Saxon potential at the appropriate binding 
energy. This potential is adjusted to generate an optimum description of the shape of the experimental cross section. 
The resulting theoretical representation of the cross section must then be multiplied by a constant factor to coincide 
with the experimental cross section. This constant factor is then interpreted as the spectroscopic factor. Another 
important ingredient in this analysis is the choice of the electron-proton cross section which must be considered 
off-shell This leads to a small additional uncertainty in the analysis of low data as discussed recently in 

Ref. [0. 

Further clarification of this intriguing situation with different spectroscopic factors at different is urgently 
needed. For this purpose we consider in this paper a study of recently published (e,e'p) data for ^^O at — 0.8 
(GeV/c)^ [|l9|. An unfactorized ap pro ach is used, as it is required at low and has been recently advocated 



also for high reactions in Ref. |20|. However, the higher energy of the outgoing proton requires a description 
which contains different elements than the conventional low analysis. For the electromagnetic current operator 
we follow the approach of Ref. where a relativistic current operator was used in a Schrodinger-based calculation, 
avoiding any nonrelativistic reduction and including the effect of spinor distortion by the Dirac scalar and vector 
potentials (see also Refs. ||2^- p5| , p^ ) . As for the final state, we employ a recently developed eikonal description of 
the final-state interaction (FSI) of the proton with the nucleus that has been tested against DWIA solutions of a 
complex spin-dependent optical potential p6| , p7t . The absorption of the proton is described theoretically by linking 
it to the corresponding absorption of a nucleon propagating through nuclear matter. The relevant quantity is the 
nucleon self-energy which is obtained from a self-consistent calculation of nucleon spectral functions including the 
effects of realistic short-range and tensor correlations |^ . This description of FSI is combined with previous results 
for the bound-state wave functions of the p-shell quasihole states in ^^O |^^, that have been deduced by solving the 
Dyson equation with a nucleon self-energy containing the same short-range correlations but in a finite volume. In 
Ref. these wave functions have been used to analyze low data for the ^®0(e,e'p) reaction In the present 
work, these very same wave functions produce a good description of the shape of the coincidence cross section for the 
p-shell quasihole states at high using the same spectroscopic factors obtained from the low data [Q. A 

preliminary result was reported in Ref. |p2| . 

A consistent analysis of low and high Q^data requires an approach in which the same quasihole wave functions and 
corresponding spectroscopic factors are used in both cases. In this paper we present such an approach. In Sec. Owe 
discuss the ingredients of the theoretical description including the structure of th e electromagnetic current (Sec. p A|), 



the eikonal approximation (Sec. II B), and the final hadronic tensor (Sec. II C). The many-particle ingredients are 
discussed in Sec. Ill which includes a summary of the calculation of the quasihole wave functions in Sec. Ill A , a 
description of the construction of the nonlocality factor required by a treatment of relativistic effects in Sec IIIB, 
and, finally , an overview of the ingredients to describe the damping of high-momentum protons in nuclear matter 
(Sec. IIIC). The results are discussed in Sec. IV while final conclusions are drawn in Sec. 



II. THE MODEL 

In the one-photon exchange approximation, the differential cross section for the scattering of an ultrarelativistic 
electron with initial (final) momentum Peipi), off a nuclear target from which a nucleon is ejected with final momentum 
p'^, reads 

^ J: l,^„w,^,,, (1) 



dPedPN 167r2 QVePe \,\'=Q,±1 



where = — uP' and q = Pe ^ p'e-, ^ = Pe~ P'e ^^^^ ^^e momentum and energy transferred to the target nucleus, 
respectively. The lepton tensor L\^x and hadron tensor W\^\i are conveniently expressed in the basis of unit vectors 



eo = (1,0,0,0) , e±i = (^0,TVi-V5*'0j ' (2) 

which define the longitudinal (0) and transverse (±1) components of the nuclear response with respect to the polar- 
ization of the exchanged virtual photon. The hadron tensor is defined as |lG,0,n3| 
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W^A,A' = ele\r E £ J,J: 5 {E, - Ef) , (3) 

I If 

i.e. it involves the average over initial states and the sum over the final undetected states (compatible with energy- 
momentum conservation) of bilinear products of the scattering amplitude . 

This basic ingredient of the calculation is built from the matrix element of the nuclear charge-current density 
operator J^^ between the initial, l^"^), and the final, \'^f), nuclear states. This complicated A-body problem can be 
simplified by projecting out of the Hilbert space the specific channel that corresponds to the experimental asymptotic 
conditions of a knocked-out nucleon with momentum p'j^ and of a residual nucleus, recoiling with momentum —pm = 
q— Ppf and mass Mr, in a well-defined state \'^^~^{En)) with energy Er and quantum numbers n. The scattering 
amplitude can be rewritten in a one-body representation (in momentum space and omitting spin degrees of freedom 
for simplicity) as [^3|jTl| ] 

Jii{u;,q,p'^,En) = I dp dp' x!,7^4„(p") J^sip, p' , q, ^) 4>e^u{p) [Sn{ER)]^ , (4) 

provided that is substituted by an appropriate effective one-body charge-current density operator J^, which 
guarantees the orthogonality between l^"^) and l^E'j) besides taking into account effects due to truncation of the 
Hilbert space. Usually, the orthogonality defect is negligible in standard kinematics for (e,e'p) reactions 10 in 
any case, J^^j is here approximated by a one-body relativistic current operator including spinor distortion along the 



lines described in the following Sec. p A 
The functions 



[Sn{ER)]-^^EnniP) - {K-\En)\a{p)\^^) , 



x[-\.M = {^t\En)\a{p)\^j) , (5) 



describe the overlap between the residual state \'^^ ^{Eji)) and the hole produced in ) and I'^f), respectively, 
by removing a particle with momentum p. Both fpERn ^^'^ Xp' "^Eru eigenfunctions of a Feshbach-like nonlocal 
energy-dependent Hamiltonian referred to the residual nucleus, belonging to the eigenvalues Er and Er + w, respec- 
tively The norm of (j^Ertn 1 SniEn) is the spectroscopic factor associated with the removal process, 
i.e. it is the probability that the residual nucleus can indeed be conceived as the target nucleus with a hole. The 
dependence of Xp' ^Eru upon p'j^ is hidden in the asymptotic state |^^) and the boundary conditions are those of an 
incoming wave. 

Because of the complexity of the eigenvalue problem in the continuum, a complex mean-field interaction with energy- 
dependent parameters is usually assumed between the residual nucleus and the emitted nucleon. Then, Xpi '^Eru x\>' ^ 
and the nonlocality of the original Feshbach Hamiltonian is taken into account by multiplying the scattering wave 
by th e appr opriate Percy factor [p5[. Several models for this FSI are discussed in the literature ( for a review, see 



Refs. |11|,|12|]). Here, the eikonal approximation is adopted and is described in more detail in Sec. [IB. Finally, in 



Sec. [I C we present the complete formula for the hadronic tensor used in the calculation. 



A. Current operator and spinor distortion 

While new data for the (e,e'p) reaction have become available at very high proton energies [^,^,|l9|, it has also 
become evident that many ingredients of the theoretical calculations must be upgraded and made adequate for the 
new kinematical regime. In particular, a nonrelativistic reduction of the electromagnetic current operator is no longer 
reliable. Since all other ingredients entering the scattering amplitude will be deduced in a Schrodinger-like framework, 
we follow the approach of Ref . |2l| . 

It is well known that a four-component Dirac spinor with positive- and negative-energy components '0+ and ?/;_, 
respectively, and satisfying a Dirac equation with energy eigenvalue E, mass m, scalar and vector potentials S and 
V , respectively, can be written as 



a-Ti D^/^{r) <j) = A(^,r) <P , (6) 

E + m + S{r) -V{r) J 
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namely it can be represented by the action of the operator A(7r, r) on the wave function (j), that satisfies a Schrodinger- 
equivalent equation with central Uc and spin-orbit Uls potentials, which can either be expressed in terms of S and 
V or replaced by intrinsically nonrelativistic potentials. The Darwin nonlocality factor 

S -V 

D(.r)^l + ^-— (7) 



is related to Uls by 



E 



ULsir) = - — — — 
2/i rD dr ' 



(8) 



where fi ~ m(A — I) /A is the reduced mass with A the mass number. 

In Ref. pH l, the calculations were performed in configuration space by defining the effective current operator 
(omitting spin indices for simplicity) 



(9) 



and by evaluating the operator part ct • vr in the effective momentum approximation (EMA), i.e. by replacing the 
operator 7? with the momenta p^,Pm determined by asymptotic kinematics. After choosing one out of the three 
(on-shell) equivalent expressions for the electromagnetic vertex function F'' the effective current operator J^jj 
was reduced to a simple 2x2 matrix acting on the nucleon spins, using the standard representation for 7 matrices as 
4x4 operators in terms of 2x2 Pauli spin matrices 

Here, the scattering amplitude is worked out in momentum space. Therefore, the operator a ■ tt becomes just a 
multiplicative factor with a 2x2 matrix structure in spin space acting on the nucleon spins. Consequently, Eq. (^) 
can be specialized to a "relativized Schrodinger framework" by considering the following effective current operator 
(omitting again spin indices for simplicity). 



J^siP,p', q, ^ / dr e*(P+^~-^^')-^ K\p' , r) 7° A(p, 



(10) 



as can be easily shown by starting from the expression of the scattering amplitude in configuration space and applying 
the proper Fourier transformations. The nonrelativistic limit of Eq. ( [lO| ) is recovered by setting E' ~ to, and 
Sir) — V{r) — 0. Inspection of Eqs. (0) and ( |6|) in dicates that in this limit the spinor-distortion operator no longer 
depends on r and the Fourier transform in Eq. (|1C|) produces the well known 5{p' — p — q) accounting for momentum 
conservation ||l^,|ri[. In the following, we will keep D{r) = 1 for the scattering state, because the distortion of a 
high-energy ejectile will be approximated by a uniform damping in nuclear matter with Uls = (see Sec. [IB). 

The electromagnetic vertex function F'^ for an on-shell nucleon can be represented through three equivalent expres- 
sions related by the Gordon identity |l^. Here, we choose the following 



pfi 

r^=7^GM(0')-:^F2(Q2) , 
2to 

where Gm is the nucleon magnetic form factor, F2 is its Pauli form factor, and 
Eq. ( pi] ) in Eq. (|lo[), the scattering amplitude becomes 



(11) 

(£" + E,p' + p). By inserting 



Jni^^q^PN^ER) 



dp dp Xp'^EEU^P) JcsiP^P'^^^) <t>ERn{P) [SniER)]' 



'dp dp' x'-^X.{p')\j 



D2 



2m 



a ■ p 



E' 



1 



2m 
a ■ p 



E' + mE- 



(27r)f 



TO 



til 



a ■ p 



E 



D- 



(7 ■ p 

E' + r. 



PA* 



a ■ p a ■ p 

E' +mE + m 



^Enn{P) [Sn{En)]-^ , 



where 



5±5 - 



{2n)i 



(12) 



(13) 



are functions of \p + q— p'\. 

The nucleon form factors are taken from Ref. while the Coulomb gauge is adopted to restore current conser- 
vation at the one-body level by modifying the longitudinal component accordingly. 
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B. The eikonal approximation 



Similar to the case of current operators, the high proton energies, that can be reached in (e,e'p) reactions at 
the new experimental facilities, also demand a suitable approach to the treatment of the proton scattering wave. 
Traditionally, the assumed mean-field interaction between the ejectile and the residual nucleus has been described by 
complex spin-dependent optical potentials with energy-dependent parameters constrained by fitting phase shifts and 
analyzing powers of elastic (inelastic) (p,p) scatterings on the corresponding residual nucleus. A Schrodinger equation 
with incoming wave boundary conditions for each partial wave of x'' , is solved up to a maximum angular momentum 

Pn 

imax(Pjv) satisfying a convergency criterion. Typically, this method has been successfully applied to (e,e'p) reactions 
at proton momenta below 0.5 GeV/c and Lmax < 50 

At higher energies, the optical analysis of proton elastic scattering is improved by the relativistic description via 
Dirac phenomenology. The scattering wave is still expanded in partial waves, but each component solves the Dirac 
equation containing the scalar and vector Dirac potentials Eo|. 

An alternative, simpler but powerful, method by Glauber j41|] suggests that, when the proton is highly energetic, 
the Schrodinger equation is reduced to a first-order differential equation along the propagation axis z, 

with boundary conditions such that asymptotically x 1, i.e. corresponding to an incoming unitary flux of plane 
waves. In the pure Glauber model, U{r) is determined in a parameter-free way starting from the elementary free 
proton-nucleon scattering amplitude at the considered energy and then averaging over all possible configurations of 
the spectator nucleons. For > 1 GeV/c, the scattering amplitude is dominated by inelastic processes and U{r) 
is supposed to be mostly sensitive to its imaginary part describing the absorption Moreover, the Glauber 

model predicts that the ratio between the real and imaginary parts of U{r) equals the ratio between the real and 
imaginary parts of the average proton-nucleon forward scattering amplitude, which is expected to be anyway small 
beyond the inelastic threshold Therefore, we can safely assume U{r) iW{r). Then, the solution to Eq. ( |l^ ) 
looks like [0 

^(-)(,^) ^ ,(-^"-^+ ^ /r U{r., z') dz') ^ J..- ^ r W[r^, z') dz') 

= eOP«-^1 e^-P^-^'^ , (15) 

i.e. as a plane wave with a damping factor related to the absorption part of the residual interaction. 

The reliability of this eikonal approximation (EA) , that has a long tradition of successful results in the field of high- 
energy proton- nucleus elastic scattering |Q] , has been tested in the context of knockout reactions and in the momentum 
range of interest here (0.6 ^ q ^ 1 GeV/c) against solutions of the Schrodinger equation with nonrelativistic complex 
optical potentials up to Lmax = 120 partial waves (see also Ref. ^]). For increasing energies, the EA is supposed 
to become more and more reliable, despite the actual semirelativistic nature of the approach Moreover, for 

emitted protons with outgoing energy beyond the inelastic threshold and initially bound momentum below the Fermi 
surface {pm ^ PFcrmi *C p'jy , (?, with pFermi the target Fermi momentum: the same kinematic conditions of the E89003 
experiment at CEBAF ||l9]), it has been shown that the proton angular distribution can actually be reproduced by 
representing the scattering wave as a plane wave with an additional damping |^.^. After all, for a fast moving 
object the nuclear density can be considered roughly constant (but for a small portion on the surface) and the eikonal 
wave of Eq. (|l5|) simply corresponds to the solution of a Schrodinger equation inside homogeneous nuclear matter. In 



the next Sec. Ill C , a microscopic justification of the damping factor will be given by a detailed description of the link 
between pi and the imaginary part of the self-energy of a nuclcon moving inside nuclear matter. Here, it is sufficient 
to say that for sake of simplicity the damping vector will be kept parallel to the wave vector of the scattered particle, 



i.e. pi II p'j^. 



The EA of Eq. (^_5|) can be also formulated by saying that the scattering wave is approximated by a plane wave 
with a complex momentum Pf = p'j^ + ipi and normalized as 

Pn ' 

with R a constant vector with modulus equal to the nuclear radius. In fact, for a propagation along the z axis, the 
wave enters the nucleus at r = — i? = (0, 0,— i?) with unitary modulus and leaves it at r = i? = (0,0, i?) damped 
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by e~'^P^'^. In order to consider the Fourier transform of Eq. (p^), an extended definition of the distribution (5 of a 
complex variable is required. In the Appendix of Rcf. ||2^ , it is actually shown that such an extension is possible so 
that we can define the EA of the scattering wave in momentum space as 



xl7\p')^e-P-^' S{pf-p') , 



(17) 



where now p' is a complex vector. The extension of the matrix element of Eq. (|T^) to the complex plane in p' is 
possible if the rest of the integrand is an analytic function asymptotic ally vanishing for \p'\ — > cx) |2^. It is rather 
easy to check that, apart from the S distribution, the integrand of Eq. (|lj) meets these requirements. Therefore, the 
scattering amplitude in the EA (also with spin indices explicitly indicated) becomes 

iJns'^n{^.q,PN.ER)^e-P-^^ J dp{s'Mip,pf,q,iu)\.Sr,) ^EnniP) [Sn{ER)]^ 

p-Pi-R 



Ef +m 



2m 



dp 



E + m 
2m 



Gm{Q^ 



Ss'^s,, + {s'n\^ ■ Pf ^ ■ P\Sn) 



{Ef + m){E + m) 

, 5i \ 



2m 



-F2iQ') 



6., 



E + m 
{s'j^\a -Pf a ■p\sn) 



a \Sr, 



Ef +m 



(Ef + m){E + m) 



(18) 



where p^ = (E,p), p'j — (Ef^pf), with E — y^lp]"^ + to^, Ef — y^\pf\^ + m^, and s'^,s„ are the projections of the 
spins of the detected proton and of the residual hole with collective quantum numbers n, respectively. The Fourier 
transform of the Darwin nonlocality factor, in agreement with Eq. (p^), is function of \p+q—pf \ — \/\p 9 ~ p'^ | ^ + Pi ■ 



C. Hadronic tensor 



After summing over the undetected final states with quantum numbers n of the residual nucleus, the hadron tensor 
Wa,a' in momentum space becomes 

W^A.A' = {-f^^'^l^V ^'^^''-^ E / dpdk {s'j,\{ls)^{pjf,q,u:) \Sn) <f>E,n{P) ^Ennik) S^^iEn) 

{sn\{Jcs)lik,Pf,q,t^)\s'N} 
EE (-)^+^'e^r e-'P-^' [ dpdk{s'^\ils)^ ip,Pf,q,io) S{p,k-En) [1^)1 {k,Pf,q,Lo) |4 ) , (19) 



where 

S{p,k-ER) = E ^"(^«) ^*Enn{k) (20) 



is the hole spectral function discussed in the next Sec. III. The isospin indices have been omitted for simplicity and, 
as before, the summation over n runs over the undetected final states of the residual nucleus that are present at a 
given excitation energy Er. 

The hole spectral function can be conveniently expanded in partial waves in a sp basis as 

S{p,k]ER) = ^ E E {^\mms\jm) {l^m[m'Jjm) Sij{p,k-,ER) Yi^n^ip) \ms){m'J Yi*^,{k) . (21) 

This expansion should not be confused with the sum in Eq. (po|): each Ij term contributes to the hadronic tensor and 
can come either from a quasihole state or from above the Fermi surface, depending on the excitation energy. 
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The angular integrations in Eq. (|T9|) can be easily performed by noting that the square root of the Darwin nonlocality 
factor, Z?^2 (r) in Eq. (^3|), is not far from 1 which would yield a 5{'p' — p — q) in momentum space (see Fig. |l]in the 



same direction as pf ~ q, i.e. 



next Sec. IIIB). Therefore, because of Eq. (0^), we impose the constraint that the vector p in Eq. (uM lies in the 



\Pf - <?! ^{p'n-^^+Pi 



and similarly for k. This approximation is reliable for high values of the involved momenta, as is the case for the 
kinematics of Ref . . It is then easy to get rid of the angular integrations in Eq. dl^ ) so that the hadronic tensor, 
with explicit spin quantum numbers, takes the form 

(W^A,A').'„ = e^v e-^'''-'^ E E E {l\mms\jm) {l\m\m',\jm)Yi„,,{pT^q)Y;^,{pT^q) 

Ij m;,mj nis.m'^ 

dpp^ i dkk^ {s'^\{J^g)^{p,pf,q,uj)\ms) Sij{p,k;ER) {m'^\{J^g)l{k,pf,q,u:)\s'j^) , (23) 
Jo 

where J^g(j),pf,q,uj) is defined by inserting the approximation ( p^ ) into Eq, 



Since the missing energy of the reaction is defined as [[13 11 



E^=Lu- Tp.^ - Tr , (24) 
where Tpi^ is the kinetic energy of the detected nucleon and 

Tr = [pfn + {Mr + Er)^] - Mr - Er (25) 
is the kinetic energy of the residual nucleus, the scattering amplitude „(w, Er) of Eq. ( p^ ) can be conveniently 

made to depend on (w, q,Pm, Em)- Therefore, the differential cross section (and other related observables) for a given 
kinematics {i-^,q) and a knockout proton corresponding to a missing energy Em will be plotted as a function of the 
missing momentum pm- Older experimental (e,e'p) data at low proton energy were usually collected in the form of 
the so-called reduced cross section ||lo| , [ll|] 

n{pm,Em) = -T^^TTT -J— ' (26) 
dPedpj^ KUeN 

where K \s a. suitable kinematic factor and (JeN is the elementary (half off-shell) electron-nucleon cross section, in 
order to reduce the information contained in a five- fold differential cross section to a two- fold function oipm a-nd E^,- 
Whenever needed, theoretical results will also be presented as reduced cross sections using the CCl prescription |l^] 
for CTeAT and the corresponding extrapolation u = E'j^~E (with E = -^Z + m?) for the off-shell nucleon. The electron 
distortion is included in the EMA by replacing q with an effective (foff @1 for the acceleration by the Coulomb field 
(in the following, for sake of simplicity the eft subscript will be omitted). 



III. QUASI-HOLE AND QUASI-PARTICLE PROPERTIES 

The calculation of the (e,e'p) cross section is most easily performed by employing the quasihole wave function of 
Eq. (^. By considering Eq. (21) in a sp basis with orbital angular momentum I, total angular momentum j, and 



momentum p, we can relate the spectroscopic amplitude to the spectral function in the following way 

Si, ip, k;E)^Y. (*o I «L, I I I O SiE {E^ Et')) , (27) 

n 

where ap/j(ajj,;j) denotes the removal (addition) operator for a nucleon. The spectral function Sij{p, k; E) can be 
obtained from the imaginary part of the corresponding sp propagator gij{p,k; E). This Green's function solves the 
Dyson equation 

gi.ip, k- E) = g\f{p, k;E)+J dpi pf J dp2 pI gf^\p,pi;E) ASi,(pi,p2; E) gi,{p2, fc; E) , (28) 

where g'^^^ refers to a Hartree-Fock propagator and ASjj represents contributions to the real and imaginary parts of 
the irreducible self-energy, which go beyond the Hartree-Fock approximation of the nucleon self-energy used to derive 
(see below). A brief summary of the calculation of the self-energy and the solution of the Dyson equation are 
included below. More details can be found in Refs. raMTIl. 
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A. Quasi-hole properties 



The self-energy is constructed by a two-step approach employing the boson-exchange potential B as defined by 
Machleidt in Rcf. |^^. The treatment of short-range correlations is taken into account by solving the Bethe-Goldstone 
equation. In the first step this equation is solved in nuclear matter at a certain density with a reasonable choice for 



the starting energy. Employing a vector bracket transformation |49|, the corresponding "Hartree-Fock" self-energy 
contribution in momentum space is calculated for ^^O using harmonic oscillator wave functions for the occupied 
states with oscillator length a = 1.72 fm~^. Since this "Hartree-Fock" self-energy is obtained from nuclear matter, 
corrections need to be applied to reinstate the properties of the ^^O Fermi surface. In addition, higher-order self- 
energy contributions are included. This procedure involves the calculation of the imaginary part of the self-energy for 
two-particle one-hole (2plh) and two-hole one-particle (2hlp) intermediate states reached by the G-matrix interaction 
and calculated in ^^O. The intermediate particle states correspond to plane waves and must be orthogonalized to 
the bound sp states [Q. Pure kinetic energies are assumed for these particle states. While this assumption is not 
very realistic for the description of the coupling to low-lying states, it is quite adequate for the treatment of tensor 
and short-range correlations. From these imaginary contributions to the self-energy one can obtain the corresponding 
real parts by employing the appropriate dispersion relations. Since the "Hartree-Fock" part was calculated in terms 
of a G-matrix, it already contains the 2plh contribution mentioned above but generated in nuclear matter. The 
corresponding real part of the self-energy as calculated in nuclear matter is then subtracted to eliminate the double 
counting terms. This procedure is quite insensitive to the original choice of density and starting energy for the nuclear 
matter G-matrix ^,^]. For the determination of the p-shell quasihole wave functions only the real part of the 
self-energy is relevant. Collecting all the contributions to this self-energy one has 

Re S]z,(p, k; E) - k) + Re Sff '(p, k; E) - Re k; E) -f Re T.]!^^''{p, k; E) 

= (p, k) + Re AE,, ip,k;E) , (29) 

with obvious notation. In the last line of Eq. ( p9| ) we have included (the real part of) ASjj which was anticipated in 
Eq. (psi). This self-energy yields a complete treatment of the effect of short-range and tensor correlations for a finite 
nucleus ||5l[] . The resulting wave functions for p-shell nucleons also yield an excellent description of the shape of the 
experimental (e,e'p) cross section |Q. The self-energy in Eq. ( p9| ) does not include an adequate description of the 
coupling of the nucleon to low- lying collective excitations which strongly influence the spectroscopic factors [^2[-p4[. 
This deficiency is not important for the present paper since we are addressing the question of the reliable extraction 
of spectroscopic factors from (e,e'p) data. It is therefore of great importance that the theoretical wave functions 
generated by Eq. ( p9| ) are of equivalent quality to the empirical ones used in the analysis of the data Indeed, 
when these wave functions are used to fit the data ||3^], they yield spectroscopic factors that are essentially identical 
to the ones from the empirical analysis. 

The solution of the Dyson equation was previously obtained in a basis generated by enclosing the system in a 
spherical box p9| , |47t . For the present paper the Dyson equation has been solved directly in momentum space by 
performing the discretization for the relevant eigenvalue problem. For discrete solutions like the quasihole states in 
^^O, the Dyson equation yields the following eigenvalue equation 

„2 poo 

f-{^^,^-^\aJ,l,\^^) + dk fc^Re Y.i,{p,k-E^){^t^auA^^) = S^^'VpyK) . (30) 
Discretizing the integration in Eq. ( |30| ) yields a straightforward diagonalization problem. The resulting quasihole 



wave functions for the pi/2 and pz/2 states are used for the analysis presented in Sec. IV. We have checked that the 



present solution method yields identical results as compared to those from the "box method." 



B. Darwin nonlocality factor 

As discussed in Sec. [I A , the Darwin nonlocality factor given by Eq. (^ is required for a proper treatment of the 
current operator at high proton energies. It is clear from Eq. (|^) that this nonlocality factor is related to the spin-orbit 
pote ntial. This relation can therefore be used to derive the nonlocality factor from the nucleon self-energy discussed in 
Sec. [II A. A complication in deriving this result is that the self-energies constructed for ^^O are inherently nonlocal in 
coordinate space. This many-body nonlocality is already present when Fock terms to the self-energy are considered. 
Additional contributions are generated when higher-order terms are included as in Eq. (p9[). The terminology here 
may be confusing so it is useful to point out that the Darwin nonlocality factor refers to the nonrelativistic reduction 
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of the Dirac equation which yields a nonlocal term in the Schrodinger equation when starting from a local Dirac 
equation. 

Since Eq. requires a local spin-orbit potential we will construct local potentials from the nonlocal self-energies. 
As shown in Rcf. [Sy], it is possible to construct local potentials from the nonlocal self-energy by using the following 
expression 



Re Sj°™'(r) 



dr' r'2 Re T.ij{r,r' ; E) , 



(31) 



where the nonlocal self-energy in coordinate space is obtained from the one in momentum space by a double Fourier- 
Bessel transformation given by 



2 /"CXj />00 

Re Si, (r, r';E)^- dp dp' p'^ ji (pr) Re S,, (p, p'; E) ji (p'r' 

Jo Jo 



(32) 



As shown in Ref. a good representation to local potentials given by Eq. ( pl| ) is generated by a Woods-Saxon form. 
Following this procedure we have obtained two local Woods-Saxon potentials V^i/2(f) and Vp3/2{r) for the relevant 
p-shell quasihole states in ^^O, respectively. Since the goal of the present work is to study the possibility to extract 
spectroscopic factors at different Q^, we have adjusted these potentials slightly to generate the correct experimental 
spin-orbit splitting. In addition, we have ensured that the corresponding wave functions in momentum space have 
the maximum overlap with those of the nonlocal self-energies. These overlaps are given by 99.97% for the Pi/2 and 
99.99% for the P3/2 wave functions, respectively. It is now possible to decompose the local potentials and V^3/2 

in central and spin-orbit potentials in the following way 



For p states this implies that 



V = Uo + UlsL ■ S 



ULs{r) = \{ypi{r)-Vpx{r) 



(33) 



(34) 



This potential can then be used to construct the Darwin nonlocality factor by inverting Eq. (||) . The latter is displayed 
in Fig. |l| together with the functions D^^(r) used in Eq. (^3|). The observed small deviation from unity of the latter 
functions is the basis for the approximation introduced in Eq. (p2|) leading to the hadron tensor (125" 



C. Damping of a quasi- particle 



As discussed in Sec. |IB|, we will assume that the damping of the nucleon on its way out of the nucleus is described 
by a corresponding process taking place in nuclear matter. Since the sp momentum is conserved in nuclear matter, 
the propagation of a nucleon through nuclear matter is diagonal in the sp momentum and can be represented by 

9iv: E) = r d. + r d. /^ii^ , (35) 

where Sp and Sh (particle and hole spectral function) describe the strength distribution above and below the Fermi 
energy for a nucleon with sp momentum p. These spectral functions have recently been determined self-consistently 
by including the effects of short-range and tensor correlations in the self-ener g y p 8 | . For this purpose the effective 
interaction is represented by the equivalent of the T-matrix in the medium [ p5|j56| . The propagation of the nucleons 
determining this in- medium interaction is also described by Eq. (|35|) which includes full off-shell effects. The resulting 
interaction is employed to construct the nucleon self-energy. This self-energy is then used to solve the Dyson equation 
for nuclear matter 

g{p;E)=9^'\p;E)+g^"Hp;E) E^^^(p;£;) g{p;E) , (36) 
where the unperturbed propagator is given by 

o(0) (v E) = - Pf) djPF - P) ,3„^ 

^ ' E-p^/2m + ir] E-p'^/2m-ir] ' ^ ' 
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FIG. 1. The Darwin nonlocality factor D{r) obtained by inverting Eq. 
D~2(r) entering Eq. ([isb, respectively. 




Dashed and dot-dashed lines show D 2 (r) and 



The solution procedure for this problem involves several iteration steps which are required because the solution to 
the Dyson equation already appears in the determination of the effective interaction and the self-energy, illustrating 
the nonlinearity of this problem. The solution of Eq. (36) can be written as 



E -jP/2m-i:{p; E) 



(38) 



Taking advantage of the slow variation of the imaginary part of the self-energy as a function of p, one can expand the 
self-energy at the momentum po for which 



Zm 



(39) 



Performing the expansion in the square of this momentum and keeping both the real and imaginary parts of the self- 
energy at Po plus the first derivative of the real part, one obtains the so-called complex pole approximation (CPA) |5|] 
for the propagator which gives a very accurate representation of this quantity when transformed to coordinate space. 
For momenta above the Fermi momentum one has 



9cpa{p; E) 



pI 



■ p^ + i"f 



where 



dp^ 



(40) 



(41) 



PS/ 



and 



7= |lmS^*'^(po;£;) 



NM 



\ dp"^ 



(42) 



This form of the propagator at a fixed energy has a simple pole structure in the complex momentum plane. The 
location of the relevant pole is given by 
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FIG. 2. The quantity 7 of Eq. (^^, related to the imaginary part of the nucleon self-energy, as function of po ~ Pjvi the 
nucleon momentum. 



Kq 



(43) 



The i mag inary part of this momentum is then used to describe the damping in the eikonal approximation described 
in Sec. [IB, i.e. ImKp = pi, with pi from Eq. ([l^). To clarify the physics further one may obtain the propagator in 
coordinate space using 



gcPA{r,r';E) = 



{2n) 



dPp e 



ip-{r—r') 



gcPAip; E) 



47r 



(44) 



From this result it is clear that the damping of the nucleon propagating through nuclear matter is determined by the 
imaginary part of kq, which in turn is determined by the imaginary part of the self-energy at this energy. In Fig. ^ 
the quantity 7 of Eq. (|4^), related to the imaginary part of the self-energy, is shown as a function of po ^ Re = p'^. 
It is remarkable that at ^ 7 fm~^, i.e. at the same proton kinematics of the NE18 experiment 7 is such 
that ImKo = Pi ^ 50 MeV/c gives the proper damping necessary to describe the observed absorption. In fact, in 
the context of the pure Glauber approximation {W oc p^) one would expect a higher proportionality factor, thus 
overestimating the quenching due to FSI (see Ref. ||2^ and references therein). The outlined derivation of 7 gives 
a microscopic explanation for reducing this proportionality factor when embedding the travelling proton in nuclear 
matter. 

This absorption effect of the medium is obtained for a realistic nucleon-nucleon interaction [5^ ] . Since this interaction 
is fitted to low-energy data, it is used in the Lippmann-Schwinger equation to describe these data. As a result, the 
coupling to intermediate states at higher energy is constrained by the fit to these low-energy data. Whether the 
description of these intermediate states as nonrelativistic two-nucleon states is accurate is then less relevant. One 
may also interpret the coupling to these intermediate (nonrelativistic two-nucleon) states as a phenomenological way 
to include the coupling to inelastic channels, quark effects, etc. For this reason we expect the present microscopic 
description of nucleon absorption in the medium to be fully relevant for the (e,e'p) reactions studied in this paper. 



IV. RESULTS 

In this section we will discuss the results for the cross section of the ^^0(e,e'p) reaction leading to the ground state 

and the first | excited state of the residual ^^N nucleus. The main theoretical ingredient is the hadronic tensor of 
Eq. (^3|), which describes the electromagnetic interaction assuming a relativistic one-body current operator including 
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spinor distortion in the initial state only. This is consistent with the spin-orbit effects associated with the quasihole 
states in the residual nucleus [see Eq. (H)]. The proton scattering wave is described in the eikonal approximation 
(EA), assuming a uniform and constant damping by nuclear matter through a nucleon self-energy containing the 
same short-range correlations used to generate the properties of the quasihole in the bound state. The electron wave 
is described through the EMA, that incorporates the acceleration due to the Coulomb field. 
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FIG. 3. Cross section for tlie ^®0(e, e'p)^^N reaction at Ep = 90 MeV constant proton energy in the center-of-mass system 
in parallel kinematics [31]. Data for state have been multiplied by 20. Solid line is the result of Ref. [30] using quasihole 
states; dashed line represents the same result but replacing the quasihole spectral function with the bound state from Ref. [58] . 
All curves have been rescaled by the spectroscopic factors Z=0.644 and Z=0.537 for the pi and p| states, respectively. 

We first reconsider this reaction at low using the conventional optical potential analysis for FSI [Q . In Fig. ^ 
the data from Ref. [^l| have been collected at a constant proton energy of 90 MeV in the center-of-mass system. 
They refer to the reduced cross section, defined by Eq. (|2^), as a function of the missing momentum pm in parallel 
kinematics, i.e. for p'j^ \\ q. Therefore, the Pm distribution can be obtained by increasing the momentum transfer 
q from positive to negative values of Pm- Two transitions were considered, leading to the ground state p^ and to 
the first excited state p| at Em = 6.32 MeV of "'^^N. The data for the transition to the p^ ground state have been 
multiplied by 20. The solid lines are the result of the calculation employing the quasihole part of the spectral function 
of Eq. (21) for the p^ and p| partial waves, respectively. The normalization of the curves is adjusted to fit the data 



indicating that the intrinsic normalization of the quasihole, 0.89 for the p^ and 0.914 for the 



reduced to Z, 



Opl/2 



/ 2 , must be significantly 

0.644 and ^op3/2 = 0.537, respectively, because only t^he depletion due to short-range correlations 

3" 



has been taken into account |3C|]. Incidentally, long-range correlations spread the total | strength over three states 
in the discrete spectrum, so that the p| data account for 86% of the strength only; by rescaling the spectroscopic 
factor by this fraction we get ^op3/2 — 0.624, in close agreement with the corresponding ground state value [^. The 
dashed lines in Fig. ^ refer to the calculation including the same spectroscopic factors but replacing the quasihole 
bound state by the wave function of Ref. [Q, which is obtained by a Skyrme-Hartree-Fock method for ^^O. Both 
descriptions are in very good agreement with the data, with a slight preference for the quasihole results at negative 

Pm- 

In Fig. U the same reaction is considered in a very different kinematical regime, namely at constant {q, uj) with 
Q2 = 0.8 (GeV/c)2 ||l^. The data here refer to a five-fold differential cross section, avoiding any ambiguity in 
modelling the half off-shell elementary cross section a^p of Eq. (|2^) . Again, results for the transition to the ground 
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^^0(e,e'p)^^N reaction at = 0.8 (GeV/c)^ in perpendicular kinematics [19]. .^^^^ 
. The solid lines represent the result of the present calculation. The dashed lines are obtained 



FIG. 4. Cross section for the 
state have been multiplied by 20 
by replacing the quasihole states with the bound state wave functions of Refs 



Data for pi 



[58]. In all cases, the results have been rescaled 



by the same spectroscopic factors as in Fig. 3, namely Z=0.644 and Z=0.537 for the pi and p| states, respectively. 



state pi have been multiplied by 20. The theoretical calculations are displayed with the same notations as in Fig. 
i.e. solid lines for the results with the quasihole bound state and dashed lines by employing the wave function of 
Ref. [Q. The preference for the first choice is here more evident. In any case, it is remarkable that the calculations 
reproduce the data by using the same spectroscopic factors as in the previous kinematics, i.e. ^opi/2 = 0.644 and 
^0p3/2 = 0.537. Therefore, contrary to the findings of Ref. [^, we do not find any need for a dependence of the 
spectroscopic factors over a wide kinematical range. This outcome is particularly welcome, since by definition these 
factors describe a spectroscopic nuclear property that must be independent of the probe scale Q^. Finally, we conclude 
that the treatment of the bound state wave function is not responsible for the dependence found in Ref. ||l6| . 

In Fig. pi we show the results for the structure functions /a, A' in the same kinematics and with the same notations 
as in Fig. M. They are defined by 

/oo = Wqq 

hi - Wi, + W'_^_^ 

/oi = 2Re [W^^ - W^^.i] 

2Re [W^^.i] , (45) 

with W'y_ y the hadronic tensor in the proton center-of-mass system. It is related to Wa,A' in the lab frame by the 
transformation Wa,a' = e*"*^'^^''' Wi^ y, i.e. by a rotation around the q direction of the angle a between the lepton 
scattering plane and the plane formed by q and p'pf. The rotation affects only the interference components, so that 
the cross section becomes pO|JTT| 

dd 1 

-pTT^ = T^Til T {-^oo/oo + iii/ii + ^oi/oicosa + cos 2a} , (46) 

i.e. it becomes parametrized in terms of the different components of the nuclear response /a, A' to the virtual photon 
probe in the spherical basis. The agreement with data is still good for both transitions over the whole Pm range, 
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FIG. 5. Structure functions for the ^^0(e,e'p)^^N reaction at = 0.8 (GeV/c)^ in perpendicular kinematics [19]. Same 
notations and scaling of curves as in Fig. 4. 



except for foi, the interference between the longitudinal and transverse responses, which is known to be particularly 
sensitive to relativistic effects ||5^,^,Q . 
Correspondingly, the left-right asymmetry 



Alt = 



da{a = 0°) - da{a = 180°) 
da{a = QO) + da{a = 180°) 



(47) 



is displayed in Fig. ^ in the same conditions and with the same notations as in Fig. ^. The discrepancy previously 
noted for /oi is here amplified, particularly for the results with the bound state of Ref. ||5|] . A possible explanation 
is that only a full account of relativistic effects, specifically of spinor distortion in both bound and scattering states, 
is needed to reproduce the data ||l^ , |59| , ^ . In the present calculation this effect is included only for the bound state, 
while the Darwin nonlocality factor for the scattering state turns out to be 1 because the homogeneous damping in 
nuclear matter does not include spin-orbit contributions. 



V. CONCLUSIONS 



We have developed a model for describing the (e,e'p) reaction at high while linking it to nonrelativistic micro- 
scopic many-body ingredients like the quasihole spectral function. The goal is to critically consider the issue raised 
in Refs. |l5| , |l^ about a possible dependence of the spectroscopic factors upon Q^. 

We use an unfactorized approach where, following Ref. [pl|, a relativistic one-body electromagnetic current operator 
is adopted in a Schrodinger-based framework avoiding any nonrelativistic reduction. The effect of spinor distortion 
by the Dirac scalar and vector potentials is consistently included only for the bound state by evaluating the Darwin 
nonlocality factor through the spin-orbit potential generated by the self-energy of the quasihole spectral functions. 
The proton scattering wave is described in an eikonal approximation (tested against DWIA solutions of a complex 
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FIG. 6. Left-right asymmetry for the '^®0(e,e'p)^^N reaction at = 0.8 (GeV/c)^ in perpendicular kinematics [19]. Same 
notations as in Fig. 4. 

spin-dependent optical potential (2^). The absorption is calculated by using a spectral function for nucleons in nuclear 
matter including the same short-range and tensor correlations adopted in the calculation of the nucleon self-energy 
in a finite volume for the p-shell quasihole states of ^^O. 

In Ref. these bound state wave functions have been used to analyze the data for ^^0(e,e'p) at low |3^] 
yielding a very good description of the reduced cross sections. In the present work, we have considered the recent data 
for the same reaction at higher pof and we have performed the analysis using the same bound state wave functions 
and the same spectroscopic factors extracted from the low-Q^ analysis. The description of the data at higher is 
still very good regarding both the 5-fold differential cross section and the structure functions. Only the interference 
/oi structure function, and the related left-right asymmetry Alt, show a visible discrepancy, particularly for the 
state. A possible explanation could be related to our incomplete treatment of the relativistic effects because the spinor 
distortion of the final state is not considered. 

However, we emphasize that our consistent analysis of low- and high-Q^ data using the same microscopic many- 
body ingredients for the quasihole states and the damping of the proton scattering wave allow us to conclude that 
we do not observe any dependence of the spectroscopic factors over the considered wide range 0.02 < < 0.8 
(GeV/c)^. This outcome is most welcome, since by definition these factors describe a spectroscopic nuclear property 
that must be independent of the probe scale . Finally, since we get a very good description of the high-Q^ data 
replacing our quasihole states with the bound states of Ref. Q , we can also conclude that the quality of the wave 
functions is not responsible for the unexpected dependence of the spectroscopic factors observed in Ref. ||l^ . 
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